Optimal Prandtl number for heat transfer in rotating Rayleigh-Benard convection 
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Numerical data for the heat transfer as a function of the Prandtl (Pr) and Rossby (Ro) numbers 
in turbulent rotating Rayleigh-Benard convection are presented for Rayleigh number Ra — 10 8 . 
When Ro is fixed the heat transfer enhancement with respect to the non-rotating value shows a 
maximum as function of Pr. This maximum is due to the reduced efficiency of Ekman pumping 
when Pr becomes too small or too large. When Pr becomes small, i.e. for large thermal diffusivity, 
the heat that is carried by the vertical vortices spreads out in the middle of the cell, and Ekman 
pumping thus becomes less efficient. For higher Pr the thermal boundary layers (BLs) are thinner 
than the kinetic BLs and therefore the Ekman vortices do not reach the thermal BL. This means 
that the fluid that is sucked into the vertical vortices is colder than for lower Pr which limits the 
efficiency of the upwards heat transfer. 

PACS numbers: 47.27.te,47.32.Ef,47.20.Bp,47.27.ek 



Turbulent Rayleigh-Benard (RB) convection, i.e. the 
flow of a fluid between two parallel plates heated from 
below and cooled from above, is the paradigmatic sys- 
tem for thermally driven turbulence in a confined space. 
Considerable progress on the understanding of the flow 
dependencies has been achieved in the last two decades 
as reviewed in [H, Q . Here we study the case where the 
sample is rotated around the vertical axis at an angu- 
lar velocity Q. That system is relevant in the astro- and 
geophysical context, namely for convection in the arc- 
tic ocean [3], in Earth's outer core [4], in the interior of 
gaseous giant planets [5[, and in the outer layer of the 
Sun H. 



Previous experimental 
merical studies 
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have shown that 
the convective heat transport increases at moderate ro- 
tation rates, before it rapidly decreases for stronger rota- 
tion. This heat transport enhancement has been ascribed 
to Ekman pumping [jj E3, El EI E EE El, ~ 
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Namely, whenever a plume is formed, the converging ra- 
dial fluid motion at the base of the plume (in the Ekman 
boundary layer) starts to swirl cyclonically resulting in 
the formation of vertical vortex tubes. The rising plume 
induces stretching of the vertical vortex tube and hence 
additional vorticity is created. This leads to enhanced 
suction of hot fluid out of the local Ekman layer and 
thus increased heat transport. Corresponding phenom- 
ena happen at the upper boundary. This process strongly 
depends on Ra and Pr and the boundary conditions 
[19|. In ref. [ll| it was shown that no heat-transfer en- 
hancement is observed at small Pr < 0.7. This was ex- 
plained by the larger thermal diffusivity at lower Pr due 
to which the heat that is carried by the vertical vortices 
created by Ekman pumping spreads out in the middle of 



FIG. 1: The ratio Nu(Q)/Nu(ft = 0) as function of Ro on 
a logarithmic scale for Ra = 1 x 10 8 and different Pr. Blue 
open circles: Pr = 0.7, red open squares: Pr = 6.4, and black 
open diamonds: Pr = 20. 



the cell. This leads to a larger destabilizing temperature 
gradient in the bulk and a lower heat transfer. Further- 
more, Ekman pumping turned out to become less effi- 
cient at higher Ra numbers where a larger eddy thermal 
diffusivity limits the efficiency of Ekman pumping. 

In the present work we systematically determine the 
heat-transfer enhancement with respect to the non rotat- 
ing value as function of the Prandtl (Pr) number by using 
direct numerical simulations (DNS). We find that in cer- 
tain regimes the heat-flux enhancement can be as large 
as 30%. Even more remarkably, we observe a heretofore 
unanticipated maximum in the heat transfer enhance- 
ment as function of Pr. 



2 



For given aspect ratio r = D/L (D is the cell diam- 
eter and L its height) and given geometry, the nature 
of RB convection is determined by the Rayleigh number 
Ra = /3gAL 3 /(kv) and the Prandtl number Pr = v j k. 
Here, (3 is the thermal expansion coefficient, g the grav- 
itational acceleration, A the temperature difference be- 
tween the bottom and top, and v and k the kinematic 
viscosity and the thermal diffusivity, respectively. The 
rotation rate Q (given in rad/s) is non-dimensionalized 
in the form of the Rossby number Ro = \J pgA/L/(2Q). 
In the DNS we solved the three-dimensional Navier- 
Stokes equations within the Boussinesq approximation. 
The numerical scheme was already described in Refs. 
3 0, El E3, 0, 0- The numerical details concern- 
ing rotating RB convection are described in detail in Ref. 



Ill ll2[ . Most simulations were performed on a grid of 
129 x 257 x 257 nodes, respectively, in the radial, az- 
imuthal, and vertical directions, allowing for a sufficient 
resolution of the small scales both inside the bulk of tur- 
bulence and in the BLs (where the grid-point density 
has been enhanced) for the parameters employed here 
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[l2| . Nu is calculated in several ways as described 
in detail in Ref. [23[ and its statistical convergence has 
been controlled. A grid of 193 x 385 x 385 nodes has been 
used for the simulations at the highest Pr number and 
to verify the results for Pr = 6.4 obtained on the coarser 
grid. Furthermore, for the higher Pr number cases the 
flow is simulated for a very long time (400 dimension- 
less time units to reach the statistically stationary state 
followed by 800 time units for averaging) to assure sta- 
tistical convergence. In refs. EH we already showed 
that our DNS results agree very well with experimental 
results in this Ra number regime. 

The numerical results for Nu(Q)/Nu(0) as function of 
Ro for several Pr are shown in figure [U The figure shows 
that the heat transport enhancement caused by Ekman 
pumping can be as large as 30% for Pr = 20. However, 
no heat transport enhancement is found for Pr = 0.7 
which is due to the larger thermal diffusivity which makes 
Ekman pumping less efficient [ll[ . In figure [T] one can al- 
ready see that the heat transfer enhancement reaches a 
maximum at a certain Pr when the Ro number is fixed. 
Indeed figure [2] confirms that the heat transfer enhance- 
ment as function of Pr reaches a maximum. Its location 
depends on Ro; namely the stronger the rotation rate 
the higher the Pr number for which the maximum heat 
transfer enhancement is found. This trend is even more 
pronounced when the heat transfer Nu(Ro) itself is con- 
sidered. 

The observation that there is a Pr number for which 
the heat transfer enhancement is largest suggests that 
there should be at least two competing effects, which 
strongly depend on the Pr number, that control the 
efficiency of Ekman pumping. One wonders why Ek- 
man pumping becomes less efficient for higher Pr num- 
bers. Clearly, it must have a different origin than the re- 
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FIG. 2: The heat transfer as function of Pr on a logarithmic 
scale for Ra = 1 x 10 8 . Black, red, blue, and green indicate 
the results for Ro = oo, Ro = 1.0, Ro = 0.3, and Ro = 0.1, 
respectively. The data obtained on the 193 x 385 x 385 (open 
squares) and the 129 x 257 x 257 grid (solid circles) are in 
very good agreement. 



duced efficiency of Ekman pumping at lower Pr. Indeed, 
there is an important difference between the high and 
the low Pr number regime, namely the relation between 
the thickness of the thermal and the kinetic boundary 
layers (BL). For the low Pr number regime the kinetic 
BL is thinner than the thermal BL and therefore the 
fluid that is sucked into the vertical vortices is very hot. 
When the Pr number is too low this heat will spread out 
in the middle of the cell due to the large thermal diffu- 
sivity. For somewhat higher Pr number the fluid that 
is sucked out of the thermal BL is still sufficiently hot 
and due to the smaller thermal diffusivity the heat can 
travel very far from the plate in the vertical vortices. In 
this way Ekman pumping can increase the heat trans- 
fer for moderate Pr. In the high Pr number regime the 
kinetic BL is much thicker than the thermal BL. There- 
fore the Ekman vortices forming in the bulk do not reach 
the thermal BL and hence the temperature of the fluid 
that enters the vertical vortices is much lower. This is 
schematically shown in figure [3l 



An investigation of the temperature isosurfaces [11 ] 
revealed long vertical vortices as suggested by Ekman- 
pumping at Pr = 6.4, while these structures are much 
shorter and broadened for the Pr = 0.7 case due to the 
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FIG. 3: a) Sketch for the low Pr number regime where the 
Ekman vortices reach the thermal BL. b) Sketch for the high 
Pr number regime where the Ekman vortices do not reach 
the thermal BL. 



larger thermal diffusivity at lower Pr JJJ . In figure [H 
we show the temperature isosur faces at Ro = 0.30 for 
several Pr to identify the difference between the high 
and moderate Pr number cases. The figure reveals that 
the vertical transport of hot (cold) fluid away from the 
bottom (top) plate through the vertical vortex tubes is 
strongly reduced in the high-Pr number regime compared 
to the case with Pr = 6.4. This is illustrated in figure 
2] where the threshold for the temperature isosurfaces is 
taken constant for all cases. Indeed, a closer investigation 
shows that the vortical structures for the higher Pr are 
approximately as long as for the Pr = 6.4 case. This con- 
firms the view that the temperature of the fluid that is 
sucked into the Ekman vortices decreases with increasing 
Pr. 

To further verify the above explanation we determined 
the horizontally averaged temperature at the kinetic BL 
height, defined as two times the height where the kinetic 
dissipation rate is highest, see Refs. 23|, 13] for details. 
This height is chosen as it indicates the position where 
effects of Ekman pumping are dominant. Furthermore, 
we also determined the horizontally averaged tempera- 
ture at the edge of the radially dependent thermal BL 
thickness, i.e. A^(r). We determined X s e l (r) by finding 
the intersection point between the linear extrapolation 
of the temperature gradient at the plate with the behav- 
ior found in the bulk, see Ref. [24|. Indeed, figure [5] 
shows that the average temperature at the edge of the 
kinetic BL decreases with increasing Pr. Thus the tem- 
perature of the fluid that enters the vertical vortices de- 
creases with increasing Pr and this limits the efficiency of 
Ekman pumping at higher Pr. Furthermore, the temper- 
ature difference with respect to the bottom plate shows 
that for the non-rotating case the kinetic BL is thinner 
than the thermal BL when Pr < 1 and thicker than the 
thermal BL when Pr > 1. The figure shows that this 




FIG. 4: 3D visualization of the temperature isosurfaces in the 
cylindrical sample of T = 1 at 0.65A (red) and 0.35A (blue), 
respectively, for Pr = 0.7 (left upper plot) and Pr = 6.4 
(right upper plot), Pr = 20 (left lower plot), and Pr — 55 
(right lower plot) for Ra — 10 8 and Ro — 0.30. The snapshots 
were taken in the respective statistically stationary regimes. 



transition point shifts towards higher Pr when the rota- 
tion rate is increased. 

Another way to identify the efficiency of Ekman pump- 
ing is to look at the horizontally averaged Re z ^ rrns (di- 
mensionless root mean square velocity of the axial veloc- 
ity fluctuations) value at the height X s e l (r) [l2[. An in- 
crease in the value of Re Zyrms shows that Ekman pump- 
ing becomes important. In figure it is shown that 
Ekman pumping is strongest for moderate Pr, because 
R^z.rms nas a maximum around Pr = 10. Note that the 
strength of the Ekman pumping measured in this way re- 
flects the measured Nu number enhancement. One can 
also deduce from figure [5] that the kinetic BL becomes 
thinner with respect to the thermal BL when the rota- 
tion rate is increased, because the temperature difference 
with the plate becomes smaller. Therefore more hot fluid 
is sucked into the vertical vortices and this further sup- 
ports the Ekman pumping effect for moderate rotation 
rates where the strongest heat transfer enhancement is 
observed. Furthermore, the maximum in the heat trans- 
fer enhancement as shown in figure [T] occurs at lower Ro, 
i.e. stronger rotation, for higher Pr. The breakdown of 
Nu at low Ro is an effect of the suppression of vertical 
velocity fluctuations through the strong rotation, which 
is seen in figure [6^. Moreover, in figure [6^ it is shown 
that the vertical velocity fluctuations are suppressed at 
higher Ro, i.e. lower rotation rate, when Pr is lower. 

To summarize, we studied the Pr number dependence 
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FIG. 5: The horizontally averaged temperature at the edge 
of the kinetic BL (circles) and at the edge of the thermal 
BL (squares) for Ra = 1 x 10 8 . Left: Blue, Red, and black 
indicates the results for Pr — 0.7, Pr = 6.4 and Pr — 20, 
respectively. Right: Black, red, and blue indicate the results 
for Ro = oo, Ro = 0.3, and Ro = 0.1, respectively. The 
results for Ro = 1 (not shown) are almost identical to the 
data obtained for Ro = oo. 
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FIG. 6: The normalized horizontally averaged rms vertical 
velocities fluctuation Re z ,rms at the edge of the thermal BL 
for Ra = 1 x 10 8 . Left: Blue circles, red squares, and black 
diamonds indicate the data for Pr = 0.7, Pr — 6.4, and 
Pr = 20, respectively. Right: Black squares, red circles, and 
blue diamonds indicate the data for Ro = 1, Ro = 0.3, and 
Ro = 0.1, respectively. 



of the heat transport enhancement in rotating RB con- 
vection and showed that at a fixed Ro number there is 
a Pr number for which the heat transfer enhancement 
reaches a maximum. This is because Ekman pumping, 
which is responsible for the heat transfer enhancement, 
becomes less efficient when the Pr number becomes too 
small or too large. At small Pr numbers the efficiency is 
limited due to the large thermal diffusivity due to which 



the heat that is sucked out of the BLs rapidly spreads out 
in the bulk. At high Pr numbers the temperature of the 
fluid that is sucked into the vertical vortices near the bot- 
tom plate is much lower, because the thermal BL is much 
thinner than the kinetic BL, and therefore the efficiency 
of Ekman pumping is lower. Furthermore, the rotation 
rate for which the heat transfer enhancement reaches its 
maximum increases with increasing Pr, because the sup- 
pression of vertical velocity fluctuations only becomes im- 
portant at larger rotation rates for higher Pr. 
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